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Angular distribution of leptons from the decay of massive vector bosons 
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We examine the transverse momentum Q± dependence of the helicity structure functions for 
massive vector bosons of mass Q in hadron reactions. We demonstrate that large logarithmic terms 
of the form ln(Q/Qx) in the helicity structure functions have the same origin as the logarithmic 
terms in the angular-integrated cross section and that they can be resummed to all orders in the 
strong coupling a s , in the same way as the angular-integrated cross section. As a consequence of 
current conservation, the resummed part of the helicity structure functions preserves the Lam- Tung 
relation as a function of Q± to all orders in a s . 



Introduction Massive vector boson production in 
hadronic interactions is a potent probe of short-distance 
dynamics in strong and electroweak interactions. The 
production of massive virtual photons, the Drell-Yan 
process tests predictions of perturbative quantum 
chromodynamics (QCD) at high orders in perturbation 
theory and constrains parton distributions. Production 
of the intermediate weak vector bosons, the W and Z, 
supplies accurate measurements of the masses of these 
bosons, and, particularly in the case of the W boson, 
provides electroweak observables that bound the mass of 
the Higgs boson within the framework of the standard 
model of particle physics Q . 

The massive virtual photon, the W boson, and the 
Z boson all have important decay modes into pairs of 
leptons. The angular distribution of these leptons, mea- 
sured in the rest-frame of the parent states, determines 
the alignment (polarization) of the massive vector bo- 
son and, consequently, more precise information on the 
production dynamics than is accessible from the spin- 
averaged rate alone. An understanding of the changes 
expected in the angular distribution as a function of the 
transverse momentum Q± is a topic of considerable im- 
portance, both for refined tests of QCD and for reduced 
systematic uncertainties on the determination of the W 
boson mass In this Letter , we derive predictions 

of the angular distribution as a function of Q±, including 
resummation to all orders in the strong coupling strength 
a s of the effects of initial state soft gluon radiation. 

The general formalism for the description of the angu- 
lar distribution in terms of helicity structure functions is 
developed for the Drell-Yan process in Ref. [B| . The cross 
section may be expressed as 
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where q denotes the four-vector momentum of the mas- 
sive virtual photon. The variables 8 and <fi are the polar 
and azimuthal angles in the virtual photon rest frame, 



and £1 is the solid angle f2(0, </>). The four independent he- 
licity structure functions Wi depend on the virtual pho- 
ton's mass Q, transverse momentum Q±, and rapidity y, 
as well as on the center-of-mass energy y/S of the pro- 
duction process. The dependence of the four Wi on Q± 
is our central focus in this Letter. 

At the most basic level, a massive virtual photon arises 
from quark-antiquark annihilation q + q — > 7* in a colli- 
sion of hadrons, and it is produced with Q± = within 
the framework of collinear factorization in perturbative 
QCD Q. First-order gluon radiation supplies finite Q±, 
through the quark-antiquark and quark-gluon subpro- 
cesses, q + q —+ 7* + g and q + g — ► 7* + q. These 
finite-order subprocesses yield singular differential cross 
sections as a function of Q± in the limit Q± — ► 0. For 
the angular-integrated cross section, da/d i q, it is well 
established that this divergence can be removed after re- 
summation of the singular terms from initial-state gluon 
emission to all orders in a s @, In this Letter, we ex- 
amine in detail the structure of the singular contributions 
for each helicity structure function Wi. We show explic- 
itly that the Q± — » singular contributions to Wt, Wt,, 
and Waa have the same origin as that of the angular 
integrated cross section and that they can be resummed 
to all orders in a s . 

Asymptotic tensor In this section, we extract the 
perturbatively singular expressions for each term in the 
Drell-Yan hadronic tensor, and we introduce an asymp- 
totic tensor that is defined to include all the singular 
terms and, in addition, to conserve the electromagnetic 
current. 

The full Drell-Yan hadronic tensor can be written as 
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W" u = - (gV - T^T") (Wt + Waa) 

~2X^X v Waa + Z»Z V {W L -W T - Waa) 
-{X^Z" +X V Z> 1 ) W A , (2) 

where = q^/Q, and X», Y M , and Z^ are orthogo- 
nal unit vectors in the virtual photon's rest frame, with 
T 2 = 1, and X 2 = Y 2 = Z 2 = -I. Different choices of 
the axes lead to different q = frames Q . In the Collins- 
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Soper frame [§], all perturbatively calculated helicity 
structure functions at order a s and beyond are singular 
as Q±/Q — > 0: and VFa have the power divergences, 
Q 2 /Ql and Q/Qj_, respectively, as well as ln(Q/Q±) 
divergences, while Wl and Waa show only ln(Q/Qj_) 
divergences [Io|, [HI- By expanding the full Drell-Yan 
hadronic tensor in Eq. ([2]) in the limit Q±/Q — ► 0, we 
obtain the following form for its singular terms 
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with two asymptotically divergent scalar functions: 
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In Eq. ©, = 6^+, = <^-, and = c^ 1 - are unit 
vectors that specify the light-cone coordinates of the lab- 
oratory frame, with n 2 = n 2 = 0, n\ = —1, n ■ n — 1, 
and n± ■ n = nj_ ■ n = 0. In this frame, the incoming 
hadron momenta are Pf — \J S/2n 11 and P^ = \/S/2 n M 



^n" andPf 
respectively, and the virtual photon momentum: 



with Q+ = y/(Q 2 +Ql)/2e* 



V(Q 2 + Qi)/2< 



and Q~ 

The singular tensor as Q±/Q — ► in Eq. |3j) is not 
current conserving since q^ ^ 0. In order to resum 

the singular terms of the hadronic tensor to all orders 
in a s , we require a tensor that includes all the singular 
terms and also conserves the current perturbatively at 
any order of a s . We use the term asymptotic tensor for 
this current-conserving tensor. We define it to be 
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It is equal to the singular tensor in Eq. ([3]) plus a minimal 
non-singular term such that = 0. 

The angular-integrated Drell-Yan cross section is ob- 
tained from the trace, do~ / d q cx —g^v W^ u . The trace of 
the asymptotic tensor in Eq. (j4|) should therefore be fixed 
by the asymptotic term W Asym in the angular-integrated 
Drell-Yan transverse momentum distribution 0]. This 



statement allows us to fix uniquely the asymptotically di- 
vergent function W^ sym in Eq. g]): JU 2 Asym = W Asym /2. 
However, the angular-integrated Drell-Yan cross section 
cannot fix the second scalar function W Asym in Eq. ((4|. 
This second function represents the singular perturbative 
behavior of the structure function We can set this 
question aside for the purposes of the present Letter and 
concentrate on the resummation of the singular terms of 
the other three helicity structure functions, Wt, Wl, and 
W AA . 

We reexpress the asymptotic tensor in terms of the 
previously defined unit vectors in the Collins-Soper frame 
as 
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Upon comparison with Eq.(2), we immediately derive the 
corresponding asymptotic helicity structure functions, 
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Our first result is this derivation that the asymptotic part 
of these three helicity structure functions is proportional 
directly to the asymptotic part of the angular-integrated 
Drell-Yan transverse momentum distribution. 

Perturbatively finite tensor We show in this sec- 
tion that the three asymptotic helicity structure func- 
tions derived in last section include all the divergent 
terms of the corresponding perturbatively calculated he- 
licity structure functions, and therefore, that we can de- 
fine a perturbatively finite tensor 
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at any order of a s . This perturbative finite tensor con- 
serves the current since the asymptotic tensor conserves 
the current. 

In terms of QCD collinear factorization [6| , we can ex- 
press the helicity structure functions in terms of parton- 
level helicity structure functions, Wi\ 



E 

a b 
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with i = T, L, AA, and incoming parton distributions 
0/ (£) °f h avor / an, i momentum fraction £. We have the 
same factorization relation for W Asym 
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For the spin-averaged 2^2 scattering processes qq 
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with x\ = Qj\fSe v and x 2 = Q/VSe v . The parton- 



to-parton splitting functions are 

1 + z 2 



Pqq{z) = C F 
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P qg {z) = T R [z 2 + (1 - z) 2 } , 
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(10) 
(11) 



with C F = {N 2 - 1)/(2JV C ) = 4/3, Tr = 1/2. Using 
Eq. ([6]) we find that as Q±/Q — > 0, the partonic-level 
asymptotic terms in Eq. @ remove all divergent con- 
tributions of the corresponding perturbatively calculated 
helicity structure functions. For the quark-antiquark an- 
nihilation subprocess, 
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For the quark-gluon subprocess 
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With Z\ and z 2 switched, Eq. (fl"3")) is also true for the 
gluon-quark subprocess. Other than the non-logarithmic 
finite piece (as Qj_/Q — ► 0) in the quark-gluon contribu- 
tion to Wl and Waa, the asymptotic tensor completely 
removes the leading term of the perturbatively calculated 
helicity structure functions as Q±/Q — > 0. 

We note that the Lam- Tung relation, Wl — 2WaAj is 
obeyed by the singular terms of the perturbatively cal- 
culated helicity structure functions. 



Polarized initial partons To gain insight into the un- 
canceled finite terms in Wl and Waa for the qg subpro- 
cess, we calculate the parton-level asymptotic and pertur- 
batively finite helicity structure functions for "polarized" 
incoming parton states, defined as the states with incom- 
ing parton polarization projected onto the difference of 

131 ]. For the quark-antiquark 
f^ ym , whereas for 



the parton helicity states 
subprocess, we find that . 
the quark-gluon subprocess 
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with 
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We note that Aw^ m is different from w^ ym in Eq. ©. 
The finite contributions in the perturbatively computed 
parton-level helicity structure functions are 
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x [AP qg (-z 2 ) - AP w (z 2 )] + 0(Qi). (16) 

The uncanceled finite term in the helicity structure 
functions Wl and Waa is proportional to 



P qg (-z 2 ) - P qg (z 2 ) = ^z 2 T R , 
for unpolarized initial partonic states, and to 



AP qg (-z 2 ) - AP qg (z 2 ) 



-4z 2 T R , 



(17) 



(18) 



for the "polarized" initial partonic states. Therefore, 
the uncanceled finite term in Wl and Waa in the limit 
Q±/Q - > can also be removed by the asymptotic term 
if the incoming quark and gluon have the same helicity 
(both positive or negative). We also observe that, at this 
order, the uncanceled term in the quark-gluon subprocess 
is proportional to the helicity flipping splitting function, 



P q - g+ (z) = P q+g -(z) =T R (l-z) 2 



(19) 



The finite term as Q±/Q — > for the quark-antiquark 
subprocess at this order is removed completely by the 
asymptotic term since the helicity flipping splitting func- 
tion for the quark vanishes at this order, P q - q +{z) = 
P q+q -(z)=0. 

Our second result in this Letter is that transverse mo- 
mentum dependent factorization for the full hadronic 
tensor, which is the basis for the Collins-Soper-Sterman 
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6-space resummation, holds at leading power in the 
Qx/Q expansion. Breaking of factorization occurs at 
subleading power but only in the helicity flipping chan- 
nel and does not supply leading logarithmic terms. 

Resummation As shown above, as Qx/Q —> 0, all 
leading large logarithmic terms in the helicity structure 
functions Wt, Wl, and IUaa are included in one asymp- 
totic function, W Asym , the same asymptotic function as 
in the resummed angular-integrated Drell-Yan cross sec- 
tion. Resummation of the large logarithmic terms of the 
Drell-Yan helicity structure functions can therefore be 
expressed in terms of the resummed contribution to the 
angular-integrated Drell-Yan cross section. Referring to 
Eq. ([6]), we obtain the resummed contribution to the he- 
licity structure functions in the Collins-Soper frame as 
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All depend on the same QCD resummed expression 
VK Resum that pertains to the angular-integrated Drell- 
Yan cross section [H, [l3[ . We remark that resummation 
deals with the logarithmically divergent terms in the per- 
turbatively calculated structure functions. There are re- 
maining finite terms as defined in Eq. |(7J) that make up 
the full contribution to the structure functions. 

Equation (|2U)l shows that the Lam- Tung relation, 
Wl = 2Waa, is satisfied by the resummed contribution 
to all orders in a s . This result is a direct consequence of 
current conservation that led us to the asymptotic tensor 
in Eq. {5]). Furthermore, since the resummed contribu- 
tion dominates the cross section in the region of small 
and intermediate values of Qx, we expect that violation 
of the Lam- Tung relation as a function of Q± should be 
relatively small. This third result of our investigation 
is consistent with what is found in perturbative calcula- 
tions at order a 2 s [14j, but it is shown here to be true to 



all orders in a s . Data also show reasonable agreement 
with the Lam- Tung relation for moderate values of Q± 

Summary and Discussion In this Letter, we inves- 
tigate the transverse momentum Qx dependence of the 
angular distribution of leptons from the decay of a mas- 
sive virtual photon produced in hadron collisions. We 
work in terms of the four independent helicity structure 
functions that characterize this angular distribution. We 
show for the first time that the Q± dependence of the he- 
licity structure functions Wt, Wl, and Waa is specified 
by the same all-orders resummed Qx distribution that 
characterizes the angular integrated cross section, up to 
Qx dependent kinematic factors that we list explicitly 
in Eq. (|2"U]) . One consequence of this work is the demon- 



stration that the Lam- Tung relation, Wl = 2U/aa, be- 
tween the longitudinal and the double-spin-flip structure 
functions is obeyed by the resummed cross sections as a 
function of Qx and that this relation is a direct result of 
the requirement of current conservation. 

In this Letter, we work entirely in the collinear factor- 
ization approach and resum large logarithmic perturba- 
tive contributions when Qx/Q — ► 0. When Q± is very 
small, one might study the angular distribution in terms 
of transverse momentum dependent parton distributions 
0. 

In further work, we intend to examine the Qx depen- 
dence of W and Z boson production, where parity vio- 
lating terms introduce additional helicity structure func- 
tions. 
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